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Побудова контуру об’ектив 1нтерактивним методом 
на основ! псевдодистаншйних карт 
нашвтонового зображення 


В этой статье предлагаются алгоритмы интерактивного выделения контура объектов на полутоновых 
изображениях на основе псевдодистанционной карты РОМ, упрощенной карты расстояний полутонового 
изображения, которая использует только те признаки изображения и объектов, которые необходимы 
для построения оптимального контура. Алгоритмы работают достаточно быстро на больших полутоновых 
снимках и позволяют получить контур высокого качества по дистанционной карте градиента полутонового 
изображения, стремящийся к линии скелета «размытой» границы объектов. 
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Введение 


Многие исследователи концентрируют свои усилия на автоматизации выделения 
объектов на изображениях с последующей ручной коррекцией полученных результа- 
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тов [1], [2]. На самом деле, полностью автоматические подходы, реализованные в ГИС, 
медицинских и инженерных приложениях зачастую дают неудовлетворительные де- 
фектные результаты, время коррекции которых вполне сравнимо с «ручным» анализом 
изображений. Таким образом, интерактивные подходы распознавания по-прежнему 
актуальны в своем развитии с тем, чтобы облегчить работу оператора при полном 
пользовательском контроле над процессом сегментации и распознавания. Следует отме- 
тить, что любой, существующий на сегодняшний день, подход, качественно решает 
задачу только для определённого класса изображений, поэтому автоматизированный 
анализ нового класса изображений диктует необходимость разработки новых или мо- 
дернизации существующих алгоритмов сегментации. 

Известно, что существуют методы интерактивной сегментации, основанные либо 
на выделении областей, либо на выделении границ [3], [4]. Инструмент «волшебная 
палочка» («таз1с \уапа юоЪ») [5], реализованный в большинстве коммерческих систем 
обработки изображений, принадлежит к первому типу и позволяет пользователю интер- 
активно выбирать затравочные (седловые) точки для дальнейшего «роста» областей, 
добавляя смежные соседние пиксели в соответствии с некоторым критерием однород- 
ности [6]. 

Метод «активных контуров» (известный как «зпаКе») — достаточно популярный 
подход для выделения границ, позволяющий устанавливать начальную кривую как 
первоначальную аппроксимацию выделяемого контура, изменяемую внешними и внутрен- 
ними «силами» на следующих этапах [1]. Окончательная граница оценивается путем 
минимизации энергии функционала. Однако этим процессом достаточно тяжело 
управлять интерактивно. Более того, реализация «таз1с \/ап4» ап4 «зпаКе» довольно 
сложна с точки зрения поддержки интерактивного визуального отклика во время вы- 
полнения алгоритма. 

Другой, широко известный подход, основанный на выделении краев, «Пуе-\лге» [3], 
для нахождения границы выделяемого объекта применяет поиск глобального опти- 
мального пути. Пиксели в алгоритме рассматриваются как узлы ориентированного 
графа, стоимость же основана на граничных признаках, приписанных дугам графа. 
Далее рассчитывается путь с минимальной стоимостью, от седловой точки к каждому 
пикселю изображения. Пользователь достигает желаемого результата сегментации, 
перемещая курсор вдоль границы объекта. Подходы такого плана, однако, могут вы- 
делять лишь площадные объекты, в то время, как множество задач требует определение 
средней линии протяженных объектов. 

В данной работе авторы предлагают два подхода для интерактивного выделения 
объектов: подход, основанный на аккумуляции связных компонент, и подход, в основе 
которого лежит парадигма «Пуе-\уте». В отличие от всех схожих технологий интерак- 
тивной сегментации, предлагаемый подход позволяет выделять как площадные объекты, 
определяя их контур, так и линейные, выделяя их скелетное представление. 

Средние линии выделяются с использованием новой схемы, обобщающей псевдо- 
дистанционные преобразования (РОТ) [7] на полутоновых изображениях, анализируя 
признаковое расстояние даже при отсутствии градиента на выделяемых областях. 


Трассировка объектов 


Как известно, цель сегментации - разделить исследуемое пространство на части 
(сегменты), обладающие сходными характеристиками или функциональными особенно- 
стями, отличными от характеристик других сегментов. Большинство успешно реализо- 
ванных в коммерческих системах — «полуинтерактивные», ориентированные на выделение 
и анализ контуров, методы, требующие вмешательства пользователя для настройки 
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параметров, для конечной корректировки результата сегментации, либо для «слежения» 
за формируемым контуром в реальном режиме времени. В целом, методика, построенная 
в соответствии с последними принципами, иначе называется трассировкой (объекта). В 
нашем случае «Масс \уапа», парадигма «Нуе-\ие» относятся к такому типу сегментации. 

Рассматривая случай, когда два линейных объекта частично перекрываются 
(рис. Та), можно предположить, что целью сегментации является векторное представ- 
ление пикселей ЕЁ] и Е? (рис. Та). В этом случае определяются два узловых пикселя, 
УТ и 52 (рис. Та), и выделяются две связные компоненты (рис. 16). Если существует 
возможность получения связного пути между 51 и 52, то векторизация может быть 
проведена автоматически без дальнейшего вмешательства со стороны пользователя. 
Однако, перекрытия и пересечения встречаются достаточно часто, следовательно не- 
обходима дальнейшая обработка для поиска пути между точками (реконструкция 
объектов). Сначала для решения данной задачи применяется утоньшение для построения 
скелета каждой компоненты (отмечен тонкой белой линией на рис. 1 и 2). Далее, опре- 
деляются характерные пиксели. Следует различать три типа характерных пикселей: 

— концевые пиксели; 

— узловые пиксели; 

— связные пиксели. 

Путь между затравочными пикселями может быть найден, как только характерные 
пиксели будут помечены. Было проведено исследование соседства каждой концевой 
точки для того, чтобы отыскать места, в которых необходимо провести линейную аппрок- 
симацию. Для каждого концевого пикселя, принадлежащего одной связной компоненте, 
подсчитывается расстояние до ближайшего пикселя на скелете другой компоненты. Пара 
пикселей, имеющих кратчайшее расстояние, соединяется линией (рис. 16). Верхний пре- 
дел К находится на пути этих соединяемых линий. В любое время пользователь может 
решить, надо ли строить линию между пикселями и как это делать: выбрать наиболее 
подходящий путь, варьировать А или переключиться в ручной режим аппроксимации. Ре- 
зультаты такой обработки (трассировки) исследуемых объектов показаны на рис 26. 


ЕЛ ЕЛ 
$1 $1 
$2 
а) 


Рисунок 1 — Взаимодействие двух объектов (а); выделенные компоненты (6) 


Рисунок 2 — Построение контура линейного объекта (а); 
и результат трассировки (6) 


При вычислении параметров, описывающих взаимное влияние пикселей изобра- 
жения, дистанционное преобразование (ОТ) является одним из наиболее известных и хо- 
рошо проработанных подходов. В определении дистанционного преобразования для 
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бинарных изображений рассматривают пиксели «объекта» (предположим, равные 1) 
и пиксели «фона» (предположим, равные 0). Для исходного изображения строится 
дистанционная карта (ОМ), в которой каждому пикселю исходного изображения 
ставится в соответствие неотрицательное число, представляющее собой минимальное 
расстояние от текущего пикселя до пикселя, принадлежащего к группе «объект». 

Свойства дистанционного преобразования широко используются для сегмен- 
тации, распознавания, быстрой математической морфологии и в других задачах обра- 
ботки изображений, однако построение полной карты для полутоновых изображений 
является сложным и затратным по времени. Предлагаемое в статье утоньшение полу- 
тонового изображения функционирует на базе предлагаемого упрощенного подхода, 
что является ключевым моментом алгоритма в целом. 


Построение псевдодистанционной карты расстояний 


Для решения проблем связных со скоростью и качеством построения скелета был 
разработан новый алгоритм, базирующийся на карте псевдорасстояния (Рзеидо- От$апсе 
Мар, РОМ), которая является упрощенной картой расстояний полутонового изображения 
и использует только особенности изображения, необходимые для построения скелета. 

Строго говоря, РОМ не является дистанционной картой, но она отражает топологи- 
ческие характеристики полугона, что достаточно для получения высококачественного 
скелета. Кроме того, алгоритм утоньшения, основанный на РОМ, намного быстрее 
обычных полутоновых алгоритмов утоньшения на растровой основе [8]. Разработан- 
ный алгоритм, включает следующие шаги: 

1) построение псевдодистанционной карты; 

2) определение основных топологических особенностей; 

3) построение полутонового скелета; 

4) бинаризация скелета. 

Дистанционная карта формируется посредством определения областей изображе- 
ния, в которых пиксели помечены их расстоянием до края объекта [9]. Для бинарного 
изображения она представлена с помощью полутонового изображения, для полуто- 
нового — ее размерность существенно увеличивается [10]. Построение карты расстояний 
для полутоновых изображений — достаточно сложная задача. Это связано с тем, что 
увеличивается размерность, полутоновое изображение трехмерное, а карта расстояний 
для него является четырехмерной. Для построения скелета не нужны дистанционные 
свойства [11]. Скелет отражает топологические свойства объекта и должен проходить 
через точки хребта, то есть через: 

— точки локальных максимумов или пиковые точки (которые не имеют соседей, 
превышающих их по величине); 

— точки седла (через которые в их окрестности пролегает единственно возмож- 
ный путь между их максимальными соседями); 

— узловые точки (единственным образом соединяющие более трех соседей, равных 
или превышающих их по величине); 

Вместе с тем, скелет не должен проходить через точки локальных минимумов 
(не имеющие соседей с меньшей величиной). 

Рассмотрим полутоновое изображение как набор слоев, когда каждый более низкий 
слой включает пиксели предыдущего слоя. Стандартное полутоновое изображение в тра- 
диционной вычислительной технике имеет 256 уровней. Таким образом, если каждый 
бинарный слой увеличить на 256 значений, будет достаточно пространства для построения 
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дистанционной карты расстояний каждого слоя (рис 3). В результате, формируется 
набор из 256 слоев изображения с диапазоном значений пикселей от 0 до 255 для каж- 
дого слоя. 
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Рисунок 3 — Представление профиля изображений: исходного полутонового (а); 
полутонового преобразованного (6); псевдодистанционной карты (в) 


Для оптимизации процесса построения РОМ используется метрика Чамфера, 
которая является простой в вычислениях и позволяет получить оптимальное соотно- 
шение скорости и качества [12]. 

Для построения карты был использован модернизированный двухпроходный 
алгоритм, на основе анализа локальных минимумов. Первый проход выполняется в 
двух направлениях: сверху вниз и слева направо. Во время прохода каждый пиксель 
меняется по следующим условиям 


(р. +Л.), если (рр+Л<руи(р:+Л>Ыр)),1=0,...п 
р=3/Кр), если (р. +Х <Щр)),1=0,...п . (1) 
р, в другом случае 
где р — значение пикселя; Г(р) — значение бинарного слоя в картинке; р;— значение 
пикселей окружения; /; — значения векторов в таблице-маске для метрики; #{— индекс 
элемента в таблице-маске для метрики, и — число элементов в таблице. 
Второй проход выполняется при тех же условиях в направлениях снизу вверх, 


справа налево. В результате формируется псевдодистанционная карта, на которой отобра- 
жены основные топологические особенности изображения (рис.4). 


Рисунок 4 — Исходное полутоновое изображение (а), полутоновое изображение, 
умноженное на 256 (6), псевдодистанционная карта (в), 
и 3О-представление псевдодистанционной карты (г) 
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Таким образом, РОМ соответствуют слоям наборов карт расстояния. Каждый 
слой РОМ начинается с уровня яркости, кратного 256, и яркость каждого пикселя 
отражает расстояние до ближайшего слоя с более низким уровнем (рис. 4). 


Определение контура объекта 


Как упоминалось выше, на сегодняшний день предложено много интерактив- 
ных методов сегментации, в том числе, основанных на глобальном поиске оптималь- 
ного пути, на ориентированном графе. 

Под оптимальным путем понимается путь с наименьшей стоимостью. 

Изображение или часть изображения рассматривается как ориентированный граф, 
в котором пиксели соответствуют узлам графа, а связи с соседями соответствуют дугам 
графа. Каждой дуге графа назначается цена или весовой коэффициент, который вычи- 
сляется по ряду свойств, характеризующих принадлежность пикселя к границе объекта. 
В описываемом алгоритме предлагается для оценки этой принадлежности использо- 
вать псевдо-дистанционную карту. 

Граф С задается множеством точек или вершин х,,х,......х, (которое обозна- 
чается через Х) и множеством линий или ребер а, ,а,,....а„ (которое обозначается 
символом А), соединяющих между собой все или часть этих точек. Таким образом, 
граф С полностью задается парой {Х, А}. 

Если ребра из А ориентированы (обычно показано стрелкой), то они называют- 
ся дугами, а граф с такими ребрами — ориентированным графом. Иногда дугам графа С 
ставится в соответствие некоторое число и’, определяемое на РОМ, называемое весом, 
или длиной, или стоимостью (ценой) дуги. 

Тогда граф С называется графом со взвешенными дугами. При рассмотрении 
пути и ‚ представленного последовательностью дуг (а,,а.......а,) , за его вес (или длину, 


или стоимость) принимают число /(11) , равное сумме весов всех дуг, входящих в [1 , т. е.: 


а > | (2) 


Для изображения / с размерами НхЙ’ каждый пиксель изображения р; опре- 
деляется как вершина графа, которая связана со своими ближайшими соседями 4; 
в четырехсвязной области (рис. 5). Каждой дуге (рь 49;) определяется весовой коэффи- 


циент "м характеризующии принадлежность пикселя к границе. Значение у? Ра, лежит 


Ру? 
в интервале [Г,Н]. Для назначения весовых коэффициентов используется пиксельное 
силовое преобразование [13]. 


Рисунок 5 — Модель ориентированного графа 


Задача поиска контура между любыми двумя точками изображения формули- 
руется как поиск пути на графе с минимальной суммарной ценой. В результате, опре- 
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деление границ объектов можно рассматривать в виде задачи оптимального пути между 
начальным и конечным узлами. Таким образом, каждому ребру графа присваивается 
свой вес, а значение общего веса, накопленного по определенному пути, должно быть 
сведено к минимуму. 

Вес определяется как значение пикселя на псевдодистанционной карте. Это самое 
простое и точное определение критерия границы, так как РОМ сохраняет топологию 
изображения. Расстояние до границы определяется по росту значений пикселей на РОМ. 
Таким образом, в рамках этого критерия, задача определения границы объекта анало- 
гична задаче утоньшения. Веса у, «, Определяются по следующей формуле: 


а, =-Р, |. (3) 
С помощью метода динамического программирования [3] можно просчитать 
минимальный вес пути между начальной и конечной точками. Следует отметить, что 
метод позволяет вычислить минимальный вес между фиксированным стартовым пик- 
селем и другими пикселями одновременно. Другими словами, только определив началь- 
ный пиксель вручную, пользователь может получить минимальный вес пути к любому 
другому пикселю. 
Таким образом, для интерактивного определения контура стартовый пиксель — 
это затравочная (седловая) точка, а конечный пиксель — есть любая «свободная» точка. 
На рис. 6 показан пример работы алгоритма для изображений разных типов. 


Ежгасйоп у. 


Ипе оГоБесё 


Рисунок 6 — Результат определения контура для: а) нейрона на цитологическом 
изображении; 6) дороги на аэрокосмическом изображении 


Выводы 


Предлагаемый алгоритм сегментации имеет эффективную (в смысле времени 
выполнения) реализацию парадигмы «Нуе-уте». Для поиска глобального оптимального 
пути предлагается использовать улучшенный алгоритм с использованием псевдодистан- 
ционной карты, позволяющей находить оптимальный путь с меньшими временными 
затратами. Также в реализации алгоритма поиска оптимального пути предлагается 
использовать эффективную схему организации внутренних структур хранения графа 
и быстрого поиска. Данная схема обеспечивает сравнительно постоянное время вставки, 
удаления и поиска узла графа с минимальной стоимостью в сортированном списке. 

В отличие от всех похожих технологий интерактивной сегментации, предлагае- 
мый алгоритм позволяет выделять как площадные объекты, определяя их контур, так 
и линейные объекты, выделяя их скелетное представление. Данное свойство было 
отестировано на задачах сегментации площадных и линейных объектов на медицинских 
изображениях и аэрокосмических снимках. 
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